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Reflection and coupling processes exhibited by plane electromagnetic waves propagated 
in an inhomogeneous horizontally-stratified anisotropic ionosphere are associated with dis- 
crete transition points or with continuous coupling regions. These arise when the fourth 
order differential equations are written in first order coupled form, and many terms in these 
equations become infinite at the transition points. This procedure is rendered more precise 
by means of a special linear transformation that reformulates the equations in a new way, 
thereby exhibiting the manner in which local coupling processes are embedded in the more 
general background process of otherwise independently propagated characteristic waves. 
To exhibit the power of the matrix algebra involved, the case of an arbitrary number of 
characteristic waves is considered; moreover, Forsterling-type coupled equations are produced 
in a more generalized form than hitherto considered, and a discussion of the equations govern- 
ing continuous coupling completes the paper. 

1. Introduction 

The equations governing the propagation of plane electromagnetic waves in a horizontally- 
stratified inhomogeneous anisotropic ionosphere may be expressed in various ways, depending 
upon the particular investigation for which the equations are required. If the 02-axis of a 
right-handed set of Cartesian axes Ox, Oy, Oz is directed vertically upwards, then it often proves 
convenient to use second or fourth order differential equations in the dependent variables E x or 
E y . Heading and Whipple [1952] and Heading [1955] have employed this representation when 
analytical solutions of the differential equations are under consideration. The reason is that 
all known analytical solutions of the field equations may be expressed in terms of generalized 
hypergeometric functions, hypergeometric functions, confluent hypergeometric functions, 
Bessel functions and so on, these functions naturally satisfying second or fourth order differen- 
tial equations. Such functions provide the reflection and conversion coefficients for a particular 
anisotropic model ionosphere, hut by inspecting the functional form of the electric field through- 
out the medium no immediate physical interpretation is possible, whereby the actual process 
of wavereflectioD and wavecoupling may be explicitly exhibited. 

Sets of equations known as coupled equations have therefore been introduced, in order to 
enable approximate solutions for the electromagnetic field to be obtained, and in order to ex- 
hihit directly the important physical processes of wave reflection and wave coupling. Forster- 
ling [1942] considered a coupled system of two second order equations; these refer to an iono- 
spheric model in which the external magnetic field is oblique but in which propagation is vertical. 
Gibbons and Nertney [1951, 1952] have used these coupled equations to obtain approximate 
numerical solutions for various ionospheric models, while Budden [1952] has used them in dis- 
cussing limiting polarization. Heading [1953] and Budden and Clemmow [1957] have given 
similar coupled equations governing the case when incidence is oblique upon an ionosphere in 
which the magnetic field is vertical. 

When both the direction of propagation and of the magnetic field are oblique, four coupled 
equations are used, each of the first order. These were first introduced by Bremmer [1949], 
but more systematically using matric techniques by Clemmow and Heading [1954]. The most 
up-to-date account of these first order coupled equations is found in chapter 18 of the recent 
text by Budden [1961], but there the fundamental simplicity and beauty of the matrix formula- 
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tion is often disguised under a multitude of symbols and approximations. Moreover, certain 
important features of the matrices occurring in the coupled equations together with the possi- 
bility of nonsingular embedding do not as yet appear to have been recognized. 

The coupled equations exhibit those ranges of height in which the propagation of the 
four characteristic waves takes place almost independently. Any point at which many of the 
coefficients in the equations become singular denotes that coupling takes place there between 
the characteristic waves concerned, the other characteristic waves not embraced by the coup- 
ling still being propagated independently. 

The present investigation consists of an examination of the equations near the coupling 
points without the necessity of introducing approximations. It will be proved that the equa- 
tions governing the coupling process at a coupling point are embedded in the equations 
governing the propagation of the waves as a whole, such that all singularities in the coefficients 
of all equations are removed at the coupling point. Moreover, although ionospheric propa- 
gation is represented in terms of four first order equations, yielding four characteristic waves 
in a homogeneous or slowly-varying medium, generalization to wave processes governed by 
n first order simultaneous differential equations yielding n characteristic waves is easily 
possible using matrix notation. The sphere of usefulness of the present investigation will 
thereby be extended if n equations are considered. 

2. Transformation of an rth Order Linear Equation 

The preliminary study of a given linear differential equation rather than a set of first 
order simultaneous linear differential equations is necessary in order to define a principal 
embedded coupled matrix. 

The rth order linear equation in the dependent variable u 



u( r) =y T -. 1 (zW r - 1) +v r _ 2 (z)u {r - 2) + . . .+v 1 (z)u'+v &)u, 



(1) 



where bracketed superscripts denote differentiation with respect to z, and where the r 
coefficients v , v ly . . ., fl r _ 1 are given functions of z, may be expressed as a linear system using 
matrix notation as follows: 
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(2) 



or 



u' = T u 



(3) 



say, where u denotes the column consisting of r elements and T denotes the square rXr matrix 
shown consisting of the coefficients occurring in (1). 

When the matrix T is diagonalized by means of the transformation 



u= A f, 



the roots of the characteristic equation 



det(T-gI) = 



are required. This equation has the explicit form 
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'Vr-ltf 



-v r - 2 q: 



-Viq—v o =0, 



possessing the r characteristic roots q u q 2} . • ., q r > 

The elements of the characteristic vector corresponding to the root qj consist of the 
cofactors of any row of the matrix T— <^I; in particular, the cofactors of the last row are 1, 
<Lh #/> • • •> <7/ -1 respectively. Hence the matrix A of the transformation may be taken to 
be 



A- 



this is an alternant matrix, whose determinant is given by 
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detA= IT (£i—qj). 

i>j 



(4) 



It follows that A is nonsingular throughout any domain in the complex z-plane in which all 
the r characteristic roots are distinct. 
Equation (3) becomes 



or 



A'f+Af' = TAf, 
f'^A^TAf-A^A f. 



The matrix A l T A is now a diagonal matrix Q, whose diagonal elements consist of the r 
characteristic roots in order ; hence 



f-Qf-A^A'f. 



(5) 



If solutions of (5) are known, then U=A f are the solutions of (3), and in particular the solution 

of (1) is 



=/i+/ 2 + . . . +/,. 



(6) 



The r eqs (5) constitute a set of first order coupled equations equivalent to the original 
eq (1). The diagonal matrix Q will be called a principal uncoupled matrix, while the terms 
A -1 A' f will be called the coupling terms, though strictly speaking coupling arises only through 
the nondiagonal elements of the matrix A _1 A'. This product A"" 1 A', whose explicit elements 
are calculated in section 5, may be called the coupling matrix. It is the presence of this coupling 
matrix that renders the r eqs (5) simultaneous and not independent, but the coupling matrix 
vanishes if the coefficients in the original eq (1) are constants, owing to the presence of deriv- 
atives in every element. The r transformed equations are then completely independent, and 
their r solutions provide r characteristic waves independently propagated when a physical 
picture is employed for their interpretation. Then 

fj=Cj exp (ftz), 

where the Cj are arbitrary constants of integration. 

When the coefficients in (1) are not constants, usually some parameter such as the frequency 
or some property of the medium attains values that cause the elements of the coupling matrix 
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to be small in magnitude compared with the elements of the principal uncoupled matrix, 
provided a domain is considered in which A is nonsingular so that its reciprocal exists. Under 
these circumstances, the nondiagonal elements in the coupling matrix are neglected, yielding r 
independent equations, whose solutions represent r characteristic waves propagated independ- 
ently through the slowly varying medium. 

In this paper, equations of the form (1) are considered only when all the r characteristic 
roots become equal at a particular point z . To this end, r new functions a, /3, . . . , p are 
introduced, defined by the r linear equations 

qj==a +k j ~ 1 P+k 2U ~ v y+ . . . +k ir - 1)(j - v P , (7) 

where k=exp (2iri/r). The matrix of the coefficients in these equations is nonsingular, 
since it is an alternant whose elements consist of the rth roots of unity 1, k, P, . . . , k T ~ l 
raised to the appropriate powers. The r new functions are therefore uniquely defined. 
Now 

r 

=ra+P S k j ~ l +y 2 k^' l> + . . . +p S fc*'-™'-". 
But if p is an integer, 

[^ k v(j - 1) l(l-k p ) = l-(k p y=l-(k r ) p ==0 ) 



so if p^O, it follows that 



Hence 



S^-^o. (s) 



a=v r __ijr. 



For example, when r=2, a and f$ are merely the two terms that arise when the quadratic 
equation is solved using the standard formula, the familiar ± signs being associated with the /3. 
When r=3, p and y arise when the reduced cubic is solved by Cardan's method. 

To lay the basis for the following sections, v r -i is now chosen to vanish identically, indi- 
cating that a=0. It will be seen in the general theory, however, that no loss of generality 
occurs on account of this assumption. Coupling is then studied by postulating that the r 
characteristic roots become equal at z=z . Under these circumstances, when a is omitted 
from (7) and when the r q's are given equal non-zero values, the rank of the matrix of the coeffi- 
cients of the right-hand side is r— 1, since the r linear expressions are dependent, while the 
rank of the augmented matrix is r. The equations are therefore inconsistent. If on the other 
hand the r equal values of the q's are zero, the first r— 1 equations (7) yield P=y= . . . =p=0. 
since the matrix of the coefficients of these first r— 1 equations is nonsingular. 

Hence, at z=z , it follows that 

q\=<L2= • • • =2r=0; 

moreover, the coefficients of (1) in which # r __i = consist of sums of products of these r roots, 
so 

*>o(zo)=«<i(2o)= • • • =«V-t(zo)=0. 
598 



This is too general for useful progress to be made, although the theory to be presented is 
equally valid without any further simplifying assumptions. The special choice is made that 



leaving a differential equation 



Vl (z)=v 2 (z) = . . . =v r _ 2 (z)=0, 



d r u , s 



(9) 



When r=2, it should be noticed that no assumption is in fact made, and it is just this case 
that usually governs wave reflection and wave coupling. The r values of the characteristic 
roots are now given by 



Another reason why these simplifying assumptions are advantageous concerns the possi- 
bility of writing down approximate solutions of (9), since (9) will form what may be termed 
a comparison equation embedded in an equation of higher order. Approximate solutions of 
this latter general equation are given in terms of approximate solutions of the comparison 
equation. When r = 2, (9) can be solved approximately in terms of the Airy integral, while 
if r > 2, solutions may be obtained in terms of functions investigated by Heading in a series 
of papers [1957 a and b; 1960], though such approximations have not as yet been published. 
But if (1) is considered with v r _i = but with no further assumptions, useful approximate ana- 
lytical solutions are almost impossible to obtain near a point z=z at which r characteristic 
roots become equal. 

Under these assumptions, matrix T becomes 
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where v is now written simply as v. 

The characteristic equation q r —v=0 has the r roots 



q j= y-i v i/r 



and the alternant matrix A factorizes thus: 
'"l ... 

v llT ... 
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say, where V and K are defined by these two square matrices. 
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The matrix K being constant, it follows that 



1 



-K- 





v~ 1/r 




-II T 



'0 

v l/r v'/rv 
2v 2/r v'/rv 



TV 



"0 
1 
2 



K 



K. 



Clearly coupling only becomes large near a point at which v'/v is singular. 

The explicit elements of this product have been calculated by Heading [I960]. They are 
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(12) 
(13) 



The coupled eqs (5) then take the form 
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f. (14) 



The first square matrix on the right hand side of (14) is the principal uncoupled matrix. 
The coupling terms become singular whenever v'/v is singular; such a point is known as a 
reflection or coupling point. Mathematically speaking, such a point is called a transition 
point. The coupling terms are singular at this point, but the original equation 

u 7 =Tu 

possesses a special matrix (10), associated in this particular case with no additional coupling 
terms. An rXr matrix of the form (10) will be called a principal embedded coupled matrix. 

3. General System 

The equations governing ionospheric radio propagation may be expressed as four linear 
first order differential equations with four dependent variables, given for example by 
Clemmow and Heading [1954]. More generally, consider n dependent variables e h e 2 , . . . , 
e n satisfying the n linear first order equations 
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Iii maxtrix notation, 

e'=Te, (15) 

where the nX^n matrix T contains n 2 elements each of which is a given function of z. In 
the ionospheric case, n=4, and some of the elements of T are zero, while others are complicated 
functions of height, the explicit forms of which would obscure the essential simplicity of the 
present argument. The matrix T is therefore more complicated than that occurring in (3) or 
the principal embedded coupled matrix (10). 

If the dependent column e is now transformed to the column f by the transformation 

e-Sf, (16) 

it follows that 

f'^S^TSf-S^S'f. 

If the n characteristic roots of T are denoted by q h q 2 , . . ., q n > let s y be the column matrix 
formed from the cofactors of T— g^I taken along any suitable row. Then the transformation 
matrix 

S=(sts 2 . . . s J 

diagonalizes the matrix T; that is, S~ ! T S = Q consists of the // roots arranged in order down 
the leading diagonal. The coupled equations become 

f / = Qf-(SiS 2 . . . sJ-^SxSa . . . S n )'f 
^Qf-S^S'f, (17) 

valid at all points at which S is nonsingular. 

If desired, a change of independent variable may easily he effected. If £=£(z), then 



yielding 



MM, 

dz dr 9 



r/f -iof-S-i' /S f 



Similarly, a change of dependent variables may be made by placing 
where the n functions <j)j(z) denotes n given functions. In matrix form, let 

f=* g , 

where $ is a diagonal matrix consisting of the n functions <j>j(z) arranged in order down the 
leading diagonal. Then 

or 

=Qg-(S*)- 1 (S#)'g, 

implying that these equations would have been produced had S<I> been used for the transforma- 
tion matrix (16) rather than S. The characteristic vector s^- is merely multiplied throughout 
by 4>j. 
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Now each cofactor of T— q f \ is a polynomial of degree n—\ at most in q h the coefficients 
being independent of the suffix j. Generally, 

"polynomial 1 of degree n— 1 at most in qj 

polynomial 2 of degree n—\ at most in qj 



polynomial n of degree n—\ at most in qj 
'"the n coefficients of polynomial 1"^ f\ 
the n coefficients of polynomial 2 



to 



the n coefficients of polynomial nj [q] l 
1 



=P 



say, where P denotes the nXn matrix formed from the coefficients of the n polynomials. The 
last column of P consists of n— 1 zero elements, while the remaining elements is unity (when 
n is even), since only one polynomial can be of degree n— 1, the others being of degree n— 2. 
Then 

r \ 1 1 . . ^ 



S=(s 1 s 2 . . . s n )=P 






ft* 



PA, 



where A denotes the alternant matrix consisting of the roots and their powers arranged in order. 
It follows that 

S^S^A-T-^PA'+P'A) 

=A 1 A'+A 1 P- 1 P'A. 

The coupled equations (17) take the form 

f^Qf-A^Af-AT^PAf, 

similar to the set (5) apart from the additional coupling terms A^P^P'Af. As before, Q 
denotes the principal uncoupled matrix, A" 1 A' is the primary coupling matrix and A -1 P _1 P' A 
is the secondary coupling matrix. 

4. Integral Equation 

Let the matrix k~" 1 A f be separated into the sum of two parts, D and (A -1 A' — D) say, 
where D is an arbitrary diagonal matrix, though this may preferably consist of certain diagonal 
terms chosen from the explicit representation of A -1 A' given in the following section. Then 

r = (Q_D)f-(A- 1 A , -D+A- 1 P- 1 P , A)f, 
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where the matrix Q — D is diagonal. 
Let f be a solution of the equation 



f' = (Q-D)f, 



this equation representing in effect n independent linear equations. Let F„ denote the nXn 
diagonal matrix containing the elements of f arranged down the leading diagonal. Consider 
the integral equation 

(19) 



f=Nf -F f* F - , (A- 1 A'-D+A- 1 P- 1 P / A)iag, 

J a 



where N denotes an arbitrary constant diagonal matrix and a an arbitrary constant. 
This is the required integral equation, for differentiation yields 

r=NIS— PJ r> - 1 (A- 1 A , -D+A- 1 P- 1 P / A)fr/z-F F - 1 (A- 1 A , -D+A- 1 P- 1 P / A)f 
=Nf;-F;F - I (Nf -f)-(A- I A'-D+A- 1 P- , P'A)f 
from (19). 
But 



FJ=(Q-D)F„ 



since all matrices are diagonal, so 



f' = N(Q-D)f„-(Q-D)(Nf -f)-(A- 1 A'-D + A- , P- 1 P'A)f 
= (Q-D)f- (A-»A'-D+ A-»P -»P'A)f, 

since N commutes with (Q— D). This demonstrates the equivalence of the differential and 
integral equations. 

5. Explicit Form of the Primary Coupling Matrix 

In order to appreciate more clearly the character of the coupling terms A _1 A'f, it is neces- 
sary (o evaluate the explicit form of the product A _1 A'. 
The reciprocal of the nth order alternant matrix 



A= 



may be calculated to have the form 
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-&(;.=()) _8o_ 

(in) (g») 



Here, S ( = 1), Si, S 2 , . . ., S n -i are the elementary symmetric functions of the n quantities 
Qii Q.2> - - -i Q.n- The symbol S p (q s =0) is used to denote the algebraic expression S p when q s is 
replaced by zero. The symbol (q p ) is used to denote the product 

(&) s (2p— ffi)&— &) • • • (&— &-i)(2p— 2*+0 • • • (&— ff») = (^/^2)«-v ( 2 °) 

where 

F=det (gl— T). 

It can then be seen that 

and 

so the general element in the product A -1 A' is given by 
(A- ] A%=f; (A-%(A% 

A-l 
11 

=S (-D-^-^ai^oXA-Dflj-^/^. (21) 

In order to simplify this expression, the polynomial of degree n— 1 is formed, whose zeros 
are q h q 2 , . . ., q n but with (^ omitted. 
This is 

=S (-1)^,(^=0) 2"-"-' 

=S (-1)"-^ B _*(2 < =0)2*- 1 , (with p=n-i) 

fc=l 

yielding upon differentiation with respect to # 

(-^~)==ib (-l) n - k S n -.«(q t =0)(k-l)(f-\ (22) 

If is^j, this result transforms (21) thus: 

(A-Ao,=Arf(-^)i 

' («<)Laff\J-2i/J 

= _£ pF/dg F_ 
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Aq = Qj 



q.j—Qi fa) 



(23) 



since F vanishes when q=qj. Out of the n—1 factors in (qj), the factor (qj — qj) obviously 
cancels with the factor (qi~qj) occurring amongst the n—1 factors of (g*). 
On the other hand, when i=j, the left-hand side of (22) may be written as 



a 5 






n (g-&), 
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(A " n '-«[iV { -4.„ 

d n "1 

=25^ log n (<z-&,) 



2i 



^ i 

being the sum of n— 1 partial fractions. 

6. Classification of Transition Points 

A specular transition point (a reflection or coupling point in physical terms) exists at a point 
z=z in the complex plane at which two or more characteristic roots attain equality. 

A transition point of order two exists at a point where q\ = q>2] such a point may be denoted 
by 2=2 12 . The first and second diagonal elements each contain one partial fraction with a 
denominator {qx — q*), so this particular partial fraction (but no others down the leading 
digonal) is singular at the transition point. 

Equation (23) shows that every nondiagnonal element is singular at z l2 along rows one and 
two only. But only columns one and two have </\ and (/'., respectively in their numerators; 
these derivatives may be singular at z n - There is, however, a difference in the order of magni- 
tude of the various coupling terms near the singularity z V2 . For the coupling coefficients 
(A _1 A') jj (i=l, 2) are larger in order of magnitude when j= 1 , 2 than when ;/ = 3, 4, . . . , n. This 
implies that coupling exists mainly only between/] and / 2 , and that / 3 ,/ 4 , . . . ,/ n are relatively 
free from coupling with ^ and / 2 . In fact, the 2X2 matrix containing the first and second 
rows and the first and second columns is the distinct factor in determining the coupled solution 
near z X2 . 

W q\ and q 2 are expressed in the form 

2i = a + 0, q 2 = a—p, 

clearly fi={) at the transition point, while the first part a=)2{qi J rq2) is a more irrelevant term 
in determining the character of the coupling at the transition point. 

Similarly, transition points of order two exist whenever q a =qb) such points may be repre- 
sented by z — z ab . Altogether, %n(n — 1) sets of second order transition points exist in the z- 
plane. The solution for /„ and/,, near z ab is locally associated with a 2X2 matrix taken from 
the a- and 6th rows and from the a- and 6th columns of A -1 A'. 

The determination just considered of the important terms governing the coupling process 
is clearly unsatisfactory, since infinite terms occur throughout. The embedding transformation 
discussed in the next sections removes this uncertainty and yields a clearer picture of the 
coupling process than the matrix A" 1 A', which merely shows where the coupling occurs. 

Every nth order equation must of necessity possess these specular transition points of order 
two, since they are mere^ given by the solutions of the well-defined equations q a = qb' But 
some equations may possess transition points of higher order, not as a matter of necessity but 
exceptionally. If a point exists at which qa^qb^qc, the point z abc is called a third order transi- 
tion point. All elements in rows a, b, c of the matrix A^A' are now singular at the point z abc , 
but only the columns a, b, c of these rows prove to be of importance, so the 3X3 matrix taken 
from these three rows and columns is the determining factor for the coupling bcteen / a ,/ 6 ,/ c . 

Similar arguments may be used to define an rth (r<n) order transition point. Such a 
point exists if q a =q b = q c — . . . at z=z abc . . ., where r of the characteristic roots attain 
equality at this point. The essential coupling between / a , f b , j Cy ... is determined by the 
r\r matrix selected from the appropriate rows and columns of A^A'. 
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These r values of the roots may be expressed according to the scheme (7), but for the 
purpose of this paper it is assumed that q j =a+k j ~ 1 fi. For the usual case when r=2, no assump- 
tion is thereby made. 

It should be noted that of the remaining n—r roots, s of these may also attain equality at 
z=z a bc • • •• If their common value 7 is not equal to a at z abc . . ., then this set of s variables is 
not coupled to the previous set of r variables. 

If various characteristic roots attain equality only as 2-^00 for a certain range of arg z, the 
coupling is described as continuous coupling rather than specular. 

Finally, irregular coupling may occur at points where the matrix P becomes singular, 
or at points where one or more of the characteristic roots become singular. It is not intended 
to discuss this form of coupling here, since it does not appear to be amenable to systematic 
treatment. 

7. Embedding Transformation 

Specular embedding is defined with respect to a particular transition point z of order r<n. 
The first order coupled eqs (18) in the n variables are to be transformed in such a way that 
the r variables associated with the transition point appear in the new equations with a principal 
embedded coupled matrix (similar to (1.0)), while the remaining n—r variables not associated 
with the transition point appear with a principal uncoupled matrix (a diagonal matrix), in 
such a way that all coupling terms in all equations are nonsingular at z . 

The independent propagation of waves not affected by the coupling is thereby exhibited, 
together with the explicit form of the rth order equation satisfied by the r coupled waves. This 
coupling is seen to be independent of the remaining n—r waves. 

The general coupled equations 

f'^Qf-A^A'f-A^P-^Af 

arc assumed to possess a transition point z at which <2i = #2= • • • =3>« If qj=a J rk j ~ 1 fi, 
where k = e 27ri/T , a and /3 being functions of z such that (3=0 at z , a new column g is defined such 
that 

f= /I, exp (fadz) °\ g=Eg 

\ U Y n — r J 

say, where \ v denotes the pXp unit matrix. Then 

Eg+Eg=QEg-A A'Eg-A^P "P'AEg, 

g'=E-QEg-E "E'g-E A ^AEg-E^A P TAEg. (24) 
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say, where K r is the diagonal matrix with elements 1, k, k 2 , . . . , k r ~ l and Q w _ r is the diagonal 
matrix with elements g r+1 , . . . , q n . 
Equation (24) now becomes 



where 



/0K r 

g' = )g-(E- 1 A- 1 A , E-fE- 1 A-'P- , P'AE)g. 

V Q„_ r y 

Finally, a new column h is denned, satisfying 

g=R-h, 



R= 



(\ 1 1 

2 kW ¥p 



. . 1 

. . k'-'p 

. . P (r -»/3 2 



J) ... I„_ J 

'B 
<0 I„_ r 

say, where B is an alternant matrix similar to (11). 

Equation (25) now becomes 

//3K r \ 
R-"h+R -h'= R 1 h-(E- 1 A- 1 A , E+E- 1 A 1 P- 1 P'AE)R- 1 h, 

V Q»_,/ 
or 

//3K r \ 
h'=R R'h-RR-'-RCE'A-A'E+E 1 A 1 P- 1 PAE)R- 1 h. 

\ o Q„_7 

The original eq (15) is transformed into this form by the transformation 

e=Sf 
= PAf 
=PAEg 
=PAER-'h 



(25) 



(26) 



(27) 



(28) 



In the next sections, it will be shown that (27) possesses the required embedding properties. 

8. Embedded Equation 
The coupling matrix occurring in (27) will be considered in two stages, firstly 

RR^'+RE^A^A'ER 1 (29) 

RE A -P P'AER 1 . (30) 



and secondly 

Since R R _1 = I, it follows that 



RR-'^-R'R 1 . 

Moreover, E and R commute owing to the special forms of E and R; hence (29) becomes 
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-RR'+E'RA'AR-'E 
= -E 'RA'AR-'E-E'RA'AR-'E 

= -E- 1 (R'A- 1 +RA-")AR- I E 
= -E-'(RA- 1 )'(AR- 1 )E 
=E- I (AR- 1 )- I (AR-')'E. 



Also matrix (30) takes the form 



(AER-O-'P-'P'CAER- 1 ). 



(31) 



(32) 



At all points where coupling may be neglected (and this includes z since the next section 
will demonstrate that all coupling matrices are nonsingular at 2 ), equation (27) becomes 



h' 


/£K r 





> 


'h 




/B 

\0 l n 


J 


c 


r 

Q 




//3BK r B" 

V o 


l 





> 



B -' 




*n — t 



The elements of h are thus separated into two distinct parts. If the first r elements are 
denoted by h r and the remaining n — r elements by h n _ r , it follows that 



"n — r ^ln — r "-n—r* 



(33) 
(34) 



It may immediately be verified that 
11 1 



£BK,B 



p kp k 2 P 
P 2 k 2 p 2 k*p 2 




Sp 

o kp 

k 2 B 



1 k-'p- 1 k~ 2 P~ 2 
1 k~ 2 p~ l k~"p~ 2 



(35) 



W . . J 



Hence the elements of h r are associated with an equation containing a principal embedded 
coupled matrix, while the elements of h n _ r are associated with a principal uncoupled matrix. 
In fact, 

hf>=Fhi=vhL, (36) 

where v vanishes at z . The separation of the principal terms has therefore been achieved. 

9. Nonsingular Coupling Terms 

Every element in the coupling matrices (31) and (32) must now be shown to be non- 
singular at z . 
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The elements of A R J occurring in (31) are considered. Since 

'111. 

Qi 22 ft 



AR -' = 



</r 



i/ ft 2 




^0 ... 

it follows that the element (AR -1 ) W when l<t<re, \<j<r, is given by 

/•(ar-%^-^- 1 . . . ji-o r 1 ^ 0-' +I 

ik- l y- 1 

Ut- r + 1 )'- 1 , 
=S aJ- 1 (A- ,,+1 )'" 1 fl"' +I 



p-i 



Bis<-iCi«'" i "'(^-w%" (,, " ,) "" i) /3- ; ' +i 

/=o L/»=i J 



Now if (/— j+1) is a multiple of /*, 

r 

v=\ 

while if (/— j+1) is not a multiple of r, this sum vanishes by (8). Hence the only terms in (3 
arising- in r(AR -1 ) - are those of the form /3 NT where N is an integer. Moreover, the lowest 
value of (7— j+1) is — r+1, so AT" can only equal 0, 1, 2, ... . Hence 

where l=Nr-\rj— 1, where A^ is a non-negative integer chosen so that 

j-l<l<i— 1. 

In other words, this element is a polynomial in f3 r . 

it follows that (AR _1 )fj contains (3 only through terms such as Nr ^ Nr ~ 1 /3 f . 

Finally, since (AR- ] )" 1 =(adj AR-^/Cdet AR" 1 ) the elements of adj AR" 1 are further 
polynomials in /3 r , while 



C02247— 61 6 
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(let AR~ 1 -(det A)/(det R) 



= n (2,-ff > )/j8 1 + 2 --- + (r - 1) detK 

from (4) and (26). Now when i and j are both less than or equal to r, the relevant factors in 
the numerator are 

(q r —q r -i) . . . (g>— gi)(gr-i— gr-2) . . . (g,_i — gi) . . . (g 2 — gi) 

which is proportional to 

jS'- 1 ^" 2 . . . 2 0. 

Hence det AR _1 does not contain (3 as a factor, so the matrix AR -1 is nonsingular at z . 

In the product (AR -1 ) -1 (AR -1 )', each element contains terms of the form P Mr (P Nr )', 
where M>0, N>0 are integers. This is proportional to pw+wr-ip^ w ' lt i x M>0, N>1, the 
lowest power of fi yielding P T ~ l fi' . If (3, when expanded in terms of z — z commences with a 
term (z—z ) b , then fi T ~~ l fi' contains the factor (z— 2 ) 6r_1 - This is nonsingular at z if br— 1>0; 
that is, if b>l/r. This implies that v(z) = fi T , when expanded as a power series in terms of 
z— z , should commence with a term that is linear at least. 

Hence the elements of the primary coupling matrix (31) are nonsingular at z under these 
circumstances. The secondary coupling matrix (32) is also nonsingular at z 0j since it has been 
shown that AR _1 is nonsingular there. It is assumed of course that P is nonsingular at z 
in the definition of specular embedding. 

10. An Explicit Case 

To deal with an explicit case, consider a second order transition point associated with the 
third order differential equation 

u" ' =v 2 u" +ViU' +v u. 

Equation (3) is 



A 

dz 



The transformation 






u= 



'() 



()" 



1 



u' 



1 1 



1 



gi 9.2 2s 



yields equation (5), which, when expanded by the results of section 5, becomes 



V = 



r q x 0^ 
q 2 

L0 qj 



f- 



gi , gi 



0.2 , g 2 



g's , gs 



qi—qs gi— g2 q\— g2 g3— gi gi— gs q_2—q\ 

gl 1 gi g2 , g2 g3 , g3 

g2— gl g3— g2 g2— gl g2— g3 g2— g3 gl— g2 

gl , gi g2 , g2 g3 , g3 

k.g3— gl g2— g3 g3— g2 gl" g3 g3— gl g3— g2j 
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If z is a point at which qi = q 2 , and if q l = a + P J q 2 =a — /3, the coupling matrix has the form 



f a'+P' , a'+0' a'- 



23 



& 



2] - q, ' 2/3 2/3 

a'+P' , a'+/3' a'-/S 



?3— gi gi— g 3 2/3 



2/3 g 3 ~ #2 2j8 g 2 — g 3 22 — q* 2/3 



k 23— 2i g 2 — g 3 g 3 — g 2 2i — 2 3 23—2i 23— 22 J 

Singularities arise on account of the terms /3', 1//3, and /3'//3. If /3 2 is approximately linear 
near z , of the form z— z , /3=(z—z)%, so the three singular terms are of the forms (z—z )~*, 
(z—z )~$ and (z—z Q )~ l respectively. The /3'//3 terms are obviously of a higher order of magni- 
tude than the others; these give rise to the particular matrix 



2/3 



1 -1 0^ 
-1 1 




but other singular elements arise in rows 1 and 2 and columns 1 and 2. 
The embedding transformation (28), in which P = I, becomes 



u 


=AER- 


h 






a 


1 




1 
~2 


<h 


22 






u 


a! 



r 

23 



"exp (fadz) 





0* 

exp (fadz) 







1 



-0- 





2j 



exp {fadz) 

a exp (fadz) exp (fadz) 

LV+/3 2 ) exp (fadz) 2a exp (fadz) 

while the coupled eqs (27) become 



n 

2s 



<lV 



using results (31) and (35) 



0) 1 > 




/3 2 


h-E-'CAR-O-^AR-O'Eh 


J> 23> 




The principal matrices yield approximately 


d /hA /0 i\/*,\ 


dz \hJ v o/W 


> 


^3=g 3 A 3 . 



exhibiting the embedding property. It should be pointed out that when approximate solutions 
are written down, it may be necessary to retain some of the coupling terms. For example, in 
the last equation, an additional term in h z must be retained on the right hand side if the W.K.B. 
approximation for h 3 is to be obtained. 
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1 1 . A Special Property of the Coupling Matrices 

Under the original transformation (16), the coupling matrix is 

A-^A'+A-T-T'A, (37) 

while under the embedding transformation (28), the coupling matrix is 

E- 1 (AR- 1 )" 1 (AR- 1 )E + (AER- 1 )- 1 P- 1 P , (AER- 1 ). (38) 

It is supposed that the last coupling matrix occurs when the first r elements of f or h are coupled 
together at z . 

However, the n—rXn—r square matrix formed by crossing out the first r rows and the 
first r columns of both (37) and (38) are identical. This may be simply demonstrated by 
partitioning all the matrices between rows r and r-\-l and between columns r and r+1, and 
expanding the products. 

Hence the n—r equations (r-\-l<j<n) in the elements of f 

/H^+ cou Pl m g terms in f u ... ,/ r +coupling terms in / f+1 , . . . f n 
and the n—r equations (r+l<j<n) in the elements of h 

AJ=g^y+ coupling terms in hi, . . . , A r +coupling terms in h T +i, . . . h n 

are identical as far as the coupling coefficients multiplying / r+1 , . . . , j n and h r+u . . . , 
h n respectively, are concerned. 

Hence any coupling diagonal terms not neglected in the equation for/J- (r+l<j<?i) are 
identical with the coupling diagonal terms not neglected in the equation for h) (r+l<j<n). 

12. Forsterling-Type Coupled Equations 

As far as ionospheric radio propagation is concerned, two differential equations each of the 
second order may be coupled together. The original form given by Forsterling [1942] and that 
considered by Budden and Clemmow [1957] are specially simple, in that the characteristic 
equation for q reduces to a biquadratic. This special feature disguises the more general formu- 
lation that is possible. 

The general formulation of this new type of coupled equation is as follows. Let z be a 
transition point of order r, namely a point at which ([i = q2= z • • • =<2V- Then the n coupled 
eqs (27) (with no singularities in the coupling coefficients at z ) may be re-expressed in the form 
of n— r+1 equations in n— r+1 variables, taking the form 

h (r r ) =P T h r J r coupling terms 
hn- r = Q n _ r hn_ r +coupling terms. (39) 

The coupling terms involve only the n— r+1 variables h r and h n _ r . This result may be 
shown merely by eliminating h u h 2 , . . . , A r -i from the last n— r+1 eqs (27), the equation 
involving h' r being differentiated r— 1 times in the process. The coupling terms now involve 
derivatives of h r , but the coefficients themselves are also derivatives, which would be small in 
a slowly-varying medium. 

Similarly, if another transition point of order s (2<s<n — r) occurs such that s of the q's 
in Q n -r become equal, the number of equations may further be reduced. 

This process is explicitly exhibited for n = 4:, such that a transition point of order 2 occurs at 
z at which qi = q2, and another transition point of order 2 at z x at which g 3 = #4- 

The given equation 

e=Te 

612 



is transformed into 



by the transformation 

where 

E = 



h'= 


r 


1 





(T 


h+ coupling terms 








P 2 

























& 













sfi 








&> 














e=PAER 1 h, 






?xp (fadz) 









(T 


, R= 


r \ l 





<P 




exp 


[fadz) 







P -P 




















1 







1 



















1„ 




J) 





1> 





Similarly, if (/3=y + 5, <74=7 — 5, and if 

F=n "Y S=p (T* 

1 1 

exp (fydz) 1 1 

J) exp (fydz) J [JO 5 -^ 

then the transformation 

e^PAER TS -j 

(in which E, R _1 , F, S _1 are mutually commutative) transforms the equation into 



y=fo l o 

(3 2 
1 



where the coupling terms are uonsingular at z and Z\. 
Now let 



j + coupling terms 



»=©■ "=©■ 



40) 



and two simultaneous matrix equations are formed involving m and n; that is, eq (40) is 

'0 



multiplied first by L Q ] Q J and then by L () () J- 

m^n+Wm+Xn, 

n'=/^ 2 ° ^ m+Ym+Zn 

\0 <5 2 / 

= Am+Ym+Zn, 
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This yields 



(41) 



(42) 



say, where W, X, Y, Z are 2X2 matrices originating by partitioning the 4X4 coupling matrix 
in (40). 

Column m is now eliminated by differentiating (42) : 

n'^Am'+A'm+rm+Ym'+Zn'+Z'n 

= (A+Y)m / + (A'+Y , )m+Zn / +Z , n 

= (A+Y)(n+Wm+Xn) + (A , +Y')m+Zn'+Z , n (from 41) 

=An+(AW+YW+A'+Y')m+Zn' + (AX+Y+YX+Z')n 
=An+(AW+YW+A , +Y / )(A+Y)- 1 (n / -Zn)+Zn / + (AX+Y+YX+ZOn 

from (42). An is the only term on the right hand side whose coefficient does not involve a 
derivative. Hence A is the principal uncoupled matrix, while all the other terms are coupling 
terms. 

The coefficient of n' is (AW+YW+A'-f Y')(A-fY) _1 . A subsidiary change of variable 
will eliminate n[ from the first equation and ni from the second equation if required. 

The above process may obviously be carried out for any number of variables. 

13. Continuous Coupling 

Let the elements of the original matrix T (see eq (15)) be functions of z that tend to 
constant limits as, say, 2-^—00 along the negative real axis. Then the n characteristic roots 
also tend to limiting values as z-^>— 00 . 

If a range of real z existed such that all the elements of T were constants throughout that 
range, then all the coupling coefficients would be zero, in which case the n coupled eqs (18) 
reduce to 

f'=Qf, 

or 

implying ^-independently propagated waves. 

More generally, if a domain of z is chosen throughout which the coupling coefficients are 
small, the coupled equations reduce to 

/H^-CA-'A'+A-'P-T'A),,/,, 

where all nondiagonal coupling terms have been neglected; such approximate equations yield 
the W.K.B.-type solutions with arbitrary constants. These approximate solutions represent 
the n characteristic waves in the medium, and they are independently propagated whenever 
coupling may be neglected. In particular, only one characteristic wave may be considered, 
namely f u with j 2 —j z = . . . ==/ w =0. For this particular solution 

e-Sf, 

so 

e\\ 6 2 * . . . 'e n —Sii' S21' . . . 'S n i, 

thereby providing in this domain the varying ratios of the values of the n original dependent 
variables associated with this characteristic wave. 

The problem demanding attention concerns what happens to these ratios as z-^—co . 

The coupling coefficients may not be neglected when: 

i. coupling is specular; the previous embedding theory yields equations dealing with this 
case; 
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ii. coupling is irregular; namely, when P -1 P' or when one of the q's is singular. This 
paper is not concerned with this ease; 

iii. coupling is continuous, namely when some of the q's lend to equality as z-^— oo. 

To be definite, let q u q 2 , . . . q r have a common limit as 2^—00 . Then the denominators 
of many of the coupling coefficients considered in section 5 tend to zero. Thus continuous 
coupling between qi, q-2, • • ■ , q r is embedded in the background provided by the complete n 
equations. Previously, specular coupling has been explicitly exhibited without approximations. 
Here, however, it proves advantageous to use approximations in the equations. 

Moreover, it is assumed that the coupling coefficients do not become large 1 in the range of 
continuous coupling. Every denominator of the form (qi — qj), l<ij<r, contains a derivative 
in its numerator, and as 2-^—00 it is assumed that the derivative is such that tin 1 whole ratio 
remains finite at least. 

Above the region of continuous coupling, a definite approximate solution is considered, 
namely f 1} with / 2 =/ 3 = . . . /„=0. Within the region, e=S f, with/ r+1 =/ r+2 = . . . =/»=0 
throughout. Then 

r r r 

ei'.€ 2 : . . . :e w =2 £u/i+2 S 2j jj : . . . : 2 S n Jj 

=£ SMtt : ± SMfi) : . . . : ± S nj (fj/fi) (43) 

j=\ j=i j=\ 

ratios depending upon the r—l values j 2 li\,fzlj\, . . • jJt/Ji- 

These ratios are denoted by r 2 , r 3 , . . ., r T respectively, and their initial values above the 
region of continuous coupling are all zero. 

The coupled equations 

f'^Qf-S^S'f 
simplify to 

/;=&/*- S (S^S'hA (l<i<r) 

i=i 

in which the diagonal elements in S^S' may oftentimes be neglected since these would be small 
compared with the q ( . 

A new set of r— 1 equations may now be derived, with r 2 , r 3 , . . ., r r as the dependent vari- 
ables. The differentiation of r, (2<i<r) yields 

1 dz\fj h /1/1 

j=i L y=i J 

being r— 1 simultaneous nonlinear differential equations for r 2 , r 3 , . . ., r r . These are the appro- 
priate differential equations governing continuous coupling. Integration starts with zero 
value's, yielding r 2 , r 3 , . . ., r r as z-^ — °o , which yield in turn the required ratios (43) of the Vs 
associated with this characteristic wave. In the region of continuous coupling, the original 
ratios S n :S 2 i: . . .:S nl are converted into the final values. 

14. References 

Bremmer, II., Some remarks on the ionospheric double refraction, Philips Research Rep. 4, 189-205 (1949). 
Budden, K. G., The theory of the limiting polarization of radio waves reflected from the ionosphere, Proc. Roy. 
Soc. A, 215, 215-233 (1952). 

615 



Budden, K. G., and P. C. Clemmow, Coupled form of the differential equations governing radio propagation in 

the ionosphere. II. Proc. Camb. Phil. Soc. 53, 669-682 (1957). 
Budden, K. G., Radio waves in the ionosphere (Cambridge Press, 1961). 
Clemmow, P. C, and J. Heading, Coupled forms of the differential equations governing radio propagation in the 

ionosphere, Proc. Camb. Phil. Soc. 50, 319-333 (1954). 
Forsterling, K., Uber die Ausbreitung elektromagnetischer Wellen in einem magnetisierten Medium bei 

senkrechter Incidenz, Hochfr. Elek. 59, 110 (1942). 
Gibbons, J. J., and R. J. Nertnev, A method for obtaining the solutions of ionospherically reflected long waves, 4 

J. Geophys. Research 56, 355-371 (1951). 
Gibbons, J. J., and R. J. Nertney, Wave solutions, including coupling, of ionospherically reflected long radio 

waves for a particular E-region model, J. Geophys. Research 57, 323-338 (1952). 
Heading, J., and R. T. P. Whipple, The oblique reflection of long wireless waves from the ionosphere at 

places where the earth's magnetic field is regarded as vertical, Phil. Trans. Roy. Soc. A, 244, 469-503 (1952). 
Heading, J., Theoretical ionospheric radio propagation (Thesis, Cambridge, England, 1953). 
Heading, J., The reflection of vertically-incident long radio waves from the ionosphere when the earth's mag- 
netic field is oblique, Proc. Roy. Soc. A, 231, 414-435 (1955). 
Heading, J., The Stokes phenomenon and certain nth order differential equations, I; preliminary investiga- 4 

tion of the equations, Proc. Camb. Phil. Soc. 53, 399-418 (1957a). 
Heading, J., The Stokes phenomenon and certain nth order differential equations, II; the Stokes phenomenon, 

Proc. Camb. Phil. Soc. 53, 419-441 (1957b). 
Heading, J., The Stokes phenomenon and certain nth order differential equations, III; matrix applications, 

Proc. Camb. Phil. Soc. 56, 329-341 (1960). 

(Paper 65D6— 163) 



616 



